Electromagnetic spin polarization on the surface of topological insulator 
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We study the spin polarization of the electrons on the surface of topological insulators under a 
dc electric field or a circularly polarized light by using Keldysh Green's function formahsm. When 
a dc electric field Ex ~ lO'^V/m is applied, a spin polarization {ay) ~ 5.2 x 10~*A~^ is induced. 
Furthermore, we also find that a light illumination induces the out-of-plane component of the spin 
polarization as a result of the inverse Faraday effect. The magnitude of the spin polarization is 
proportional to the square of the lifetime r and (ctz) ~ 2 x 10~^°A~^ for typical parameters. 
Finally, we investigate the spin polarization in the presence of the warping term. By the symmetry 
consideration of the system, we find that the out-of-plane spin polarization is cubic of the current 
and that the magnitude of the induced spin polarization depends on the direction of the applied 
current. 

PACS numbers: 78.20.Ls, 72.10.-d, 78.68.-fm 
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I. INTRODUCTION 

The three-dimensional topological insulator (TI) is an 
intriguing material, which is insulating in its bulk but 
has conducting surface states topologically protected by 
the time-reversal symmetryji"— The surface state of the 
TI forms Dirac cone band, and recent spin and angular 
resolved photoemission spectroscopy (spin-ARPES) ex- 
periments have revealed the Dirac band structure. Also, 
other experimental reports on the TI include the obser- 
vation of a massive Dirac fermion in magnetically doped 
TIs.— Up to now, several materials, such as Bi2Te3 and 
Bi2Se3 have been identified as TIsi^i"— 

There are several theoretical predictions on phenom- 
ena unique to the TI.-'^'^^'^'^ Among them, it is worthy of 
attention that the spin and the momentum of an elec- 
tron are locked to each other, which causes external- 
field-driven spin-related responses on the surface of the 
TI. This effect will not only feature spintronics appli- 
cationsii"— but also bring up the possibilities of a new 
probe to surface states of TIs. From this point of view, 
it is worthwhile to focus on the response of the surface 
states of TIs to external electric fields. 

The spin-related response of the TI to electric fields 
is attributed to the spin-orbit coupling (SOC). A similar 
effect is also present in the Rashba system;2^ The Rashba 
system has the Hamiltonian of the form ^ -\- a{(T x fe)^, 
where Ci's are the Pauli's matrices in spin space and a is 
a constant unique to the material. Since the spins of the 
electrons in the Rashba system are aligned in the clock- 
wise or counterclockwise manner in the fc-space, it is seen 
that dc currents induce spin accumulations.—""— Because 
of this fact, we expect that interesting spin-related phe- 
nomena also appear on the surfaces of the TI. 

Magneto-optical Kerr and Faraday effects of a TI thin 
film have been studied in Refs. [27ll28l |. The inverse Fara- 
day effect (IFE) is an effect that circularly polarized light 
induces the stationary magnetization.—""— Recent exper- 
imental progress has made it possible to reverse mag- 
netization in magnets with a circularly polarized light 



through the IFE,^^""^^^ and now the IFE is a powerful 
method of ultrafast magnetization manipulation. In view 
of this, it is quite timely to study the IFE of the TI. 

Additionally, it has been clarified that the effect of 
the warping term may be remarkable in Bi2Tc3.'^^ The 
warping term arises accordingly to the space group of 
the crystal. Theoretically, based on the k ■ p theory it is 
found that the warping is due to the cubic term of k in 
the Hamiltonian.— 1^ It is easily seen that on the warped 
band structure, the spins of the eigenstates have out-of- 
plane components. Therefore it is naively expected that 
dc current flowing through the surface will induce out- 
of-plane spin polarizations. 

In this paper, we study the spin polarization of the 
electrons on the surface of TIs under a dc electric 
field or a circularly polarized light by using Keldysh 
Green's function formalism. When a dc electric field 
Ex ^ lO^V/m is applied, a spin polarization {oy) ~ 
5.2 X lO^^A^'^ is induced. Furthermore, we find that 
a light illumination induces the out-of-plane component 
of the spin polarization as a result of the IFE. The magni- 
tude of the spin polarization is proportional to the square 
of the lifetime r and {oz) ~ 2 x 10~^°A~^ for typical 
parameters. Finally, we also investigate the spin polar- 
ization in the presence of the warping term. By the sym- 
metry consideration of the system, we find that the out- 
of-plane spin polarization is cubic of the current and that 
the magnitude of the induced spin polarization depends 
on the direction of the applied current. 

The organization of this paper as follows. In section 
II, we explain the formalism. In section III, we present 
the calculated results of the current-induced spin polar- 
ization and the spin polarization by circularly polarized 
light, and discuss the spin polarization in the presence of 
the warping term. In section VI, we summarize our re- 
sults. In Appendices A and B, details of the calculation 
of Green's functions are presented. In Appendix C, we 
compare our results with those in the Rashba system and 
in Appendix D, calculation based on the Kubo formula 
is presented. 
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II. FORMULATION 

We consider a surface of a TI along the xy plane. The 
Hamiltonian of the surface state is given by 



Ho = v{(T X k). 







ivk- 




where k± — iky, cr^'s are the Pauli matrices, v is 

a constant representing the velocity, and we choose the 
unit so that h = 1. The eigenvalues are Eks = svk 

with s = ± and 



k — \ /l? 



1.2 



and the corresponding 
/2l2(l, ~isk+/k) with 



eigenstates are 0fc,s(^) — e /y 
I? being the size of the system. 

In what follows, we introduce the Keldysh Green's 
function formalism, which allows us to calculate physical 
quantities in non-equillibrium states. At first, we define 
the lesser Green's function with the field operators c and 

as G^{r,t,r' ,t') = i{cl^{r',t')cii{r,t)), which takes a 
2x2 matrix form since these field operators are spinors. 
The is the Heisenberg representation of the annihila- 
tion (creation) operator. The brackets (• • • ) mean both 
quantum and thermal averages. 

Usually, it is impossible to calculate the Green's func- 
tion exactly for complicated Hamiltonian, and therefore 
the pertubative method for the Green's function is em- 
ployed. We now consider the entire Hamiltonian written 
in the form H = + V, where the Hq is the unpertur- 
bative Hamiltonian, which has been solved exactly, and 
the 1^ is a perturbative term. The calculation is done 
with the Dyson's equation. With go and g being the un- 
perturbed and perturbed Green's functions, respectively, 
the Dyson's equation yields 

9^^9^+9oVg<+g<Vg\ (2) 

f^Q and are the retarded and advanced Green's func- 
tions, which are expressed as g^ = {uj — Ho + iO)^^ and 
gQ = {oj — Ho — iO)~^ respectively. We utilize rela- 
tions between these Green's functions go , go and go ■ 
9ok,uj = fi^)(9ok,oj - 9ok,J with /(w) being the Fermi 
distribution function4^ Then we can use Eq. ([2]) itera- 
tively to obtain arbitrary higher order perturbations as 



5o< 



9hV9^ 



- 9^Vg^ + 9lVglVg< + g^Vg^Vgl 



9^Vg-oVg-o + ---.^ 

Here, we consider spin-independent impurities V = 
X^il^T'' ^mp'5(T — Ti), where A^imp is the number of the 
impurities in the system. Supposing that the impurities 
are dilute, the self-consistent equation can be solved per- 
turbatively. We take the random average over all possible 
configulations of the impurities to recover the translation 
symmetry of the system. The self-consistent Born ap- 
proximation is equivalent to the Born approximation in 
the limit of dilute impurities. We assume this limit in 
the following, and then we get the self-energy as 



r(a) _ 



'imp ''imp 
~~L2 



9k,i 



^imp^mp r(a) 
■ 2^ 9ok,u> 
k 



L2 



limpl^mp 
4^2 



(3) 



where (Tq is the identity matrix in spin space and nimp = 
Ainip/L^ is the concentration of the impurities. We here 
supposed IefIt ^ 1 and the real part of the self-energy 
can be neglected. Using the self-energy, we obtain the 
modified Green's functions: g'*^^^ = (u — Ho — Y7'^^'^)~^. 



(1) We define rj^} — TiimpV^mpl'^l/^^^^; and the relaxation time 



is then given by r = l/2rj^^. In the following sections, 
we assume that [cfIt ^ 1 and fir ^ 1, where SI is the 
frequency of the external field. The explicit form of the 
modified Green's function reads 



r(a) 

9k,LJ 



{uj ± irj^y - u2fc2 \ ivk+ uj 



Lo —ivk- 



(4) 



which corresponds to the diagram in Fig. [TJ 



\q + 



X X 

-q,/ 'q, ''92 
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FIG. 1: The perturbed Green's function obtained by the self- 
consistent Born approximation. The single and double lines 
mean the unperturbed and perturbed Green's functions, re- 
spectively, and the dashed lines are scattering potentials by 
the impurities. 



III. RESULTS 



A. Current-induced spin polarization 



In this section, we will study the spin polarization in 
response to the dc electric field applied parallel to the 
surface. When we consider the dc external field, it makes 
the calculation easier to introduce the field oscillating 
with frequency and then take the limit Q \, -1-0. In 
the following, we consider the spatially uniform field so 
that the wave vector k becomes a good quantum number. 
We omit the wave number k of quantities for brevity. 

The external field is expressed by the vector poten- 
tial A{t) which satisfies E{t) = —A{t), and then the 
perturbation term reads V{t) — —j ■ A{t), where j is 
the current operator. The current operator is given by 
jx — —ex = evay, jy = —ey = —evox- When the elec- 
tric field is along the +x direction, the perturbation term 
yields V{t) = evExC^^^'ay / (i^) , whose Fourier transform 
is 27r(5(n - fl')VQ> = 2TrS{fl - n')evEx(Ty/ {iVl'). Using 
Eq. ([5]), we obtain up to linear order in V\ 
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Sin')g<+5in - {gLVn'gl - (r o a)) + 



n 



' {/(^) {gLV^riglf - (r O a)) + /'(c.)5LV^o' (ffS " 5^)} 



(5) 



r 



where /(w) = (g('^-«F)/fcBT _^ j^^-i -g ^j^g Fermi distri- 
bution function. We here expanded the Fermi distribu- 
tion function and Green's functions in terms of fl. In 
the above expansion, we can show that the gj VnigJ 
terms vanish by the integration in fc-space (See Appendix 

A). Furthermore, the gl^'^^Vnigl^'^^)'^ terms do not con- 
tribute to physical quantities. The higher order deriva- 
tives of Green's functions give higher order terms of 57 
and can be neglected. Then, only the 5^Vh'5^ term con- 
tributes to the result. Therefore, the Keldysh formalism 
within the linear response is apparently equivalent to the 
Kubo formula approach. Calculation based on the Kubo 
formula will be also discussed in the Appendix D. 

Next, we calculate the sum of the ladder diagrams in 
Fig. [21 which corresponds to the vertex correction. 



Vn should be doubled Vb = Vn + En=i = 2Vh. Af- 
ter substituting Vq with Vh and the Fourier transform, 
we obtain the perturbed lesser Green's function including 
the vertex correction at zero temperature of the form; 



G<(r,t,r,t) 



-y 

k 



2^ 



nf'{u;)glVngt 



(6) 



The local spin density of the spin polarization (ai) at 
the position r and time t is given by making the trace of 
the product of the Pauli's matrix and the lesser Green's 
function: 



(a,) = -iTT(7,G<ir, t,r, t) 



(7) 



where i — x,y^ z. Therefore, the expectation value of Oy 
reads 




ks 



k' s' 




iV (n)=\/ (n+1) 
' n O 



FIG. 2: We have to include the diagrams shown in (a) to 
satisfy the Ward-Takahashi identity, which ensures the gauge 
invariance. (b) and (c) are the iterative method we use to get 
the n-th ladder diagrams shown above. 



To calculate these ladder terms, we use the iterative 
procedure. First, we consider the diagram with impu- 
rity scattering in Fig. [2Jb). The corresponding term 

is (Vimpgi;Vh5S^mp)AV, where the (• • • )av is the ran- 
dom average over the impurity configurations. Using the 
concrete expressions of Green's functions and perturba- 



- T/(l 



' — 

y c — 



tion terms, we obtain (Knp5'L^a5S^mp)AV 
Vo/2. This calculation shows that the first-order cor- 
rection for the vertex function becomes half the orig- 
inal one. Similarly, the higher correction terms yield 

(Vimp5^V's^")gSVimp)AV ^ - ^s^"V2, therefore 



(") 



Vj2/2". Summing up all the ladder diagrams, 

= Vh. This means 
that when the ladder terms are considered, the original 



we obtain Y.°^^^ V^'^ 



ievEx 



evEx 



TrnimpVjjjjp 



(8) 



In a similar way, the expectation values of and tr^ are 
calculated to be zero. As for the charge conductivity, it 
can be calculated in a similar way: 



TmimpKmp 



(9) 



Note that the conductivity is independent of the Fermi 
energy within the self-consistent Born approximation, 
with the real part of the self-energy neglected. The form 
of the conductivity (Eq. ([9|)) is the same as that of the 
graphene for the short-ranged impurities^-- which also 
has a similar linear dispersion. This relation can be un- 
derstood in an intuitive way by noting a difference of the 
degrees of freedom: Graphene has two Dirac cones at K 
and K' points in its Brillouin zone, and each Dirac cone 
is spin-degenerate, and hence is two-fold. Therefore, the 
number of freedom of the graphene that contributes to 
the conductivity is four times greater than that of the 
TI. However, the contribution from non-magnetic impu- 
rities in graphene is four times larger than that in the 
TI because of the spin- momentum locking. As a result, 
these factors cancel out each other and we get the same 
expression of the conductivities. 
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Experimentally, the presence of the bulk conduction 
carrier would hinder precise estimations of longitudinal 
conductivity of surfaces of TIs. The bulk conduction is 
mainly due to vacancies of Se and Te atoms, ionic impu- 
rities or lattice defects. However, the spin polarization 
discussed above is a surface effect, and its measurement 
give us the physics of helical surface states separated from 
the bulk. In fact, there is a relation between the spin re- 
sponse function ipy and the longitudinal conductivity Uxx 
as a XX = sv(j}y, which offers a possibility of a new probe 
to the surface states of TIs. 

Now, we make a realistic estimate of the spin polariza- 
tion. We here suppose that it is a good approximation 
that the conductivity on surface of TI large as 

that of graphene; crf^'^P^™^ = 2.5x IQ-^f^-i. The param- 
eter V in the model Hamiltonian (Eq. ([ij) is w ~ 2.0eVA 
from the observation of surface states of Bi2Se3 by the 
ARPES experiment which is typically larger than the 
Fermi velocity in graphene. We assume applying the elec- 
tric field of strength lO^V/m, and then the induced spin 
polarization is calculated to be [ay) ~ 5.2 x 10~^A~^. 
The corresponding magnetization can be obtained by 
multiplying it by the Bohr magneton /ib and the elec- 
tron spin g-factor. 



being real, and Vq, and Vj] are reduced to be 



2ev ( 



~ t _ 2e« / £ 

n loo 



(10) 



Because we are especially interested in the dc response 
of the TI to the oscillating field, we calculate the second 
order perturbation with the Keldysh formalism, and take 
into account terms in Fig. 21 



oj-n 

(a) 



(b) 



FIG. 4: (a) and (b) correspond to the terms gujVgui-nVgui 
and guiV gui+nV Qui respectively. These terms give the diagonal 
components in the Fourier transformed cj space. 



B. Inverse Faraday effect 

In this section, we investigate the nonlinear effect 
caused by a circularly polarized light, namely the IFE, 
where a dc spin polarization is induced by a light illumi- 
nation. The axes are set as in Fig. |3l 



By using the Dyson's equation iteratively, we obtain 
the expression up to the second order perturbation with 
respect to V: 

G< ~g<+g'Vg<+g<Vg- 
+ 9'Vg'Vg< + g'^Vg<Vg- + g<Vg'^Vg\ (11) 



circularly polarized light 




topological insulator 



FIG. 3: (color online) The circularly polarized light propa- 
gating along the —z axis induces the spin polarization (ctz). 



In particular, the diagonal elements in uj are necessary 
to calculate the dc response of the system. is given 
by 



9^+f{^)[{9lVn9l+,M{gl-9l) 

+ f{u + mM9l+n - 9l+nM9l 

-f J{lo - n)givl{gl_^ - gl^^Wngl, (12) 



We employ the gauge to set the scalar potential to 
be zero and hence A{t) = ~E{t). Then, the perturba- 
tion term is V{t) = j ■ A{t) + h.c. We consider spa- 
tially uniform electric field E{t) = E^e''^^ with E^i = 
{E^\ Eq \ 0), then Fourier transform of the perturba- 
tion Hamiltonian is Vfc,^' — 2TTL^Sk,o[Sii^'—^)Vn+5{fl'+ 
ri)VQ], where Vn = j ■ E^/ifl and is the hermitian 
conjugate of Vq. When the polarization is left circularly 
polarized, the electric field is Eq = £{1, —i, 0) with £ 



where we used the formula g^ = f{'^){9^~9L)- Notably, 
the first perturbative terms, which include a single Vh 
, do not contribute to the result, because they do not 
contribute to the diagonal uj components. 

To proceed further, we assume that f2 is small and ex- 
pand the Fermi distribution function and Green's func- 
tions in terms of f2. We expect that we have to expand at 
least up to the third-order with respect to 51 because the 
perturbation terms Vq and Vq are inversely proportional 
to Q. Actually, the terms up to second order perturba- 
tion in become zero (See Appendix B). The leading 
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terms in Eq. (|12l) are therefore reduced to 



q2 a— r a2 a-r 

„r T> _yu; f>t„a „r T>t Jtj a 



r7 9gt- 



did 



dgtr\y „a 



dijj 



V^g't 



+^/"'(^) [glVngl-'V^,9l glV^,gl-'Vngl 



(13) 



where — g^—gl,- Second and higher derivative terms 



of the Fermi distribution function is partially integrated, 
which yields only the first order derivative of the Fermi 
distribution function. With the partial integration and 
the Fourier transformation, we get the perturbed Green's 
function G'^{r,t,r,t). From this Green's function, we 
calculate the local spin density as 



(a,) = -iTT[a,G<{r,t,r,t)] - '^Kn, (14) 



where kp/An is the 2D electron density, and is rep- 
resented with 7 = 77ep/eF as 



J 



Kn = 



27(3 - 7^) 27(3 + 272 + 37'') + 3(1 - 72)(1 + 72)2(7rsgn(7) - 2 arctan7) 



97r(l +7^)3 



127172(1+72)2 



(15) 



Kii is the averaged spin polarization per electron. Note 
that I7I - |(eFr)-i| < 1. In the limit of 7 ^- 0, Eq. (fTC)) 
gives 

^o=isgn(7) 3 , (16) 



which is proportional to the square of the relaxation time 
T. In this limit, the local spin density (cr^) yields 



(ctz) = sgn(eF)- 



(■Eq X E*^),T^n 



(17) 



Notice that a finite spin polarization per unit area, 
(ctz), is generated as a result of the light illumination. 
The spin polarization is proportional to {En x E^)z and 
hence the IFF is dependent on the helicity of the applied 
electric field. It is seen that when the polarization of the 
incident light is reversed or equivalently we change the 
field from En = f (1, — i, 0) to Eq — £{l, i, 0), then the 
induced spin polarization reverses. Apparently (ctz) is 
zero for a linearly polarized light. Although the IFF also 
occurs under the illumination of elliptic-polarized lights, 
it is most efficient in the case of circular polarized lights. 
For example, if the external field is Eq — £(1, —i, 0), 
then E^^ = £(1, +i, 0) yields En x £;* = 2i£2e^ with 
Bz being the unit vector along the z axis. As for the 
frequency dependence, the spin polarization by the IFF 
is proportional to il and vanishes in the limit of il \ 0. 

The polarization is proportional to the square of the 
relaxation time t. This is in stark contrast to the case 
of the Rashba system, where it is proportional to t~2 
for high light frequency.— The resulting spin polarization 
is independent of the sign oi ep- Therefore, for a left- 
circularly polarized light, we have —i^En x Ein)z > 
and (cr,) > 0. 



At the end of this section, let us estimate the spin 
polarization for the realistic parameters; 1/2t ~ ImeV, 
n ~ ImeV, ep ~ 102meV and v ~ 2eVA^^. When the 
electric field is of strength 5 ~ 1 x lO^V/m, then Kn = 
1 X 10^^. For the typical surface state of the TI, the 2D 
electron density is 2 x lO-^A then the spin polarization 
density is {(Tz) ~ 2 x 10~^°A~2. As noted in the former 
section, we can translate this value to the magnitude of 
magnetization by multiplying it by Bohr magneton /jb 
and the electron spin g-factor. 



C. Warping effect 

In this section, we investigate effects of the warping 
term. In terms of the crystal structure, Bi2Te3 has the 
rhombohedral structure^ and its band structure has Cgy 
symmetr y. T his effect have been treated theoretically in 
Refs. j4oll4l| . According to these studies, the effective 
Hamiltonian includes the hexagonal warping term and 
takes the form 



H ^ v{o- X k)z + -{kl + kl)az. 
The eigenstates and the eigenvalues are given by 



(18) 



^ikr 



^2 



^ik-r 



ie"l"' cos ^ 



ek^s = S\/ i;2fc2 + Q!2fc6 (,Qg2 30^^ 



(19) 



(20) 



with = tan ^{ky/kx) and tanflfc = fc2/(fc2 cos30fc), 
where kc — yjv/a. It is apparent that each eigenstate 
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has nonzero expectation value of <Jz as 

{ks\az\ks) — scosOk ^ 0, 
which is schematically depicted in Fig. [S] 



(21) 
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FIG. 5: (color online) The spin configuration of the eigenvec- 
tors for the s = -I- band. The arrows in this figure indicate the 
in-plane spin components of the eigenstates and the solid cir- 
cle (cross) indicates the out-of-plane spin components. At the 
cusp points the spin lies completely in the xy plane just as in 
the case without the warping term. However, at other points 
the spin has nonzero component, which is the largest in its 
magnitude at the valley points. 



To investigate the current induced spin polarization, 
we expand the spin density (ct^) in terms of the current 
(ji) = ^tU) with coefficients Aij: 



(22) 



where the expression Aij is restricted by the symmetry 
of the system. 

The Hamiltonian (Eq. (HH])) possesses the time-reversal 
symmetry T, the threefold rotational symmetry around 
the z-axis Csz, and the mirror symmetry with respect 
to the yz-plane, My^. The mirror symmetry My^ gives 
{jx) -{jx), (jy) (jy) and (ct^) -{<Jz), and there- 
fore ((Tz) is an odd function about (j^) and Aij is allowed 
to be nonzero only for i — odd. In addition, because of 
the time-reversal symmetry, (az) is an even function of 
(jy) and Aij is allowed to be nonzero only for j — even. 
Furthermore, (ctz) is invariant under the three-fold ro- 

tation Csz-. iUx), (jy)) ^ (-i(i.) - ^{jy), ^{jx) - 
i(jj,)), which gives a restriction on the coefficients Aij. 
For example, i = 1, j = term can not be invariant 
under the 27r/3 rotation, but the i = 1, j = 2 and 
« = 3, J = terms remain the same under this rotation 
if = —A12/3. From these symmetry considerations, 
the expansion of the spin polarization with respect to the 
current in the lowest order yields 



{cTz)^{jx)^~mx){jy)' 



(23) 



which is cubic of the current. When the external field 
is small, the current is linearly dependent on the field. 



and the spin polarization can be written as (ctz) oc — 
SE^Ey. The spin polarization depends not only on the 
magnitude but also the direction of the electric field. The 
spin polarization is largest when the field is along the x 
axis, and zero when the field is tilted by 7r/6 from the x 
axis. 



IV. CONCLUSION 

In this paper, we have investigated the response of the 
surface of TI to the electric fields. We have used the 
Keldysh formalism to calculate the spin polarization of 
the surface of TI in a systematic fashion. In the linear re- 
sponse regime, we have found that dc current can induce 
the spin polarization perpendicular to the current due to 
the strong spin-orbit coupling. Also, we have estimated 
the magnitude of the spin polarization. We have further 
investigated the inverse Faraday effect of the TI based on 
the second order perturbation, and we have found that 
the spin polarization emerges by circularly polarized light 
illumination, which is proportional to square of the life- 
time r. We also investigated effects of the warping term 
and found that there is a spin polarization perpendicular 
to the surface of TIs which has the three fold rotational 
symmetry and is cubic with respect to the current. 
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Appendix A: Calculation of the Green's function in 
Eq.(5) 

Here, we show how the g'^^'^^fVg'''^'^'' terms vanish in 
Eq. ([5]). The key is that V is written in the form of cur- 
rent operator. Hence, the vertex function can be written 
in the form of the derivative of the Hamiltonian in terms 
of k r'. 



ieEx dHo 
Q dkx ' 



and hence we have 



r(a)T> r(a) ieE^ dg'^"-^ 

g ^ 'Vng ^ ' = 



(Al) 



(A2) 



Considering the integration of these terms about fc, we 
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obtain 



the derivative of the Hamiltonian with respect to k± : 



1 

4^ 



kx — +00 



(A3) 



where the integrand is zero from the periodicity of the 
system. Therefore, the integration becomes zero after 
alL 



Appendix B: Calculation of the Green's function in 
Eq.(12) 



O dk+' 



2ie£ dHo 
n dk-' 



(Bl) 



Therefore, we have 



r2 aA:4 



dk_ 



(B2) 



Here, we show details of the calculation of the Green's 
function in Eq.(12). In this case, V^^^ is proportional to 



The terms independent of and composed of only g^''^^ 
in the Eq. are shown to be vanish as follows 



d^ 
4^ 



{g'Vng'V^^g' + .g'Vo/Vh/ - (r o a)) cx 



d2fc 

4^ 

d2fc 

4^ 



5/ a/ 



-1 dg^ 



dk+dk- 



-g" + 9 



.dig! 



dk- 

^ dg" dg^ 



(r o a) 



9fc_ dk. 



+ 



dk. 



(B3) 



where we used the partial integration and made use of the 
relation d{g')-^g' /dk± = 0. Note also that (gL^^'V^ ^ 
id — Hq ± irji^ and the Hamiltonian Hq is linear about k. 



Next, we show that the terms linear in which consist 
of only vanish as follows: 



d^ 



g'Vn^V^g' - gX%Vng^ - ^ ^) 



cx —f{uj)fl 



0, 



d2fc 

4^ 
fk 
47r2 



g'Vng'g'V^g' - g'V^g'g'Vng' - (r o a) 



dg'^ dg" dg^ 



i9fc+ dk- dk- dk+ 



r o a 



(B4) 



where we again used the partial integration. As for those tion function, we can take a similar procedure to show 
terms which include the derivative of the Fermi distribu- that they vanish: 
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= 0, 



d2fc 

4^ 



dk+ dk. 



-9 



^dHodg^^dgldHo , 
dk- dk+ dk- dk+ 



dk^ 



' dk+ 



9+9 



d fdHo 



dk+ \ dk^ 



-9 



^ dHo dg^ ^ dg' dHo ^ 
^ dk^ dk+ dk^ dk+ " 



47r2 



dk-dk+' 



^ dk+dk- 



(B5) 



since d^Ho/dk+dk- = 0. 

The terms proportional to or fi^ have to be zero be- 
cause otherwise some physical quantities, such as electron 
density oc TrG^(r, r, t), diverge in the low frequency 



limit n\0. 

Regarding the terms proportional to fl^, we find that 
these vanish after multiplying the corresponding term by 
az and taking trace: 



9LVn^vU' 



9lvf^Vn9l 



h.c. 



a2 a— r a-r 

9lVn-l^vU + 9lV^^:^Vn9l 



9u,^n 



duj 



n9u. 



9LV^^^Vngl 



= 0, 
(B6) 



with g**"'' = g"" - g"- 



Appendix C: Comparison to the Rashba system 

The previous studies on the Rashba systems^^ have 
shown that the dc current-induced spin accumulation oc- 
curs in the Rashba system^ir— . In this Appendix, we in- 
troduce the calculation of longitudinal conductivity and 
the spin accumulation in the Rashba system. The calcu- 
lation methods used below are presented in Ref. [l^. 

The model Hamiltonian for the Rashba system is 

^Rashba (CI) 

where the a is the magnitude of the SOC. The eigen- 
values are given by Ep-s = k^ /2m + sak, where k = 

^ k"^ + ky and s = ±, and the corresponding eigen- 

states are 0s(^) = e*'"''/\/2L^(l, —isk+/k). It is no- 
table that the Rashba Hamiltonian has the diagonal 
component in its Hamiltonian in contrast to the case 
of the TI. The current operator is — e{hkx — otUy) 
and jy = e(bky + aa^) with b = 1/m. The unper- 
turbed Green's function is gQ — [z — ^Rashbaj-i_ taking 



into account spin- independent short-range impurity po- 
tential V{r) = V^X]fc^r^^(^ ~ the self-energy can 
be derived. The modification can be made via the self- 
energy in a similar manner to the case of the TI. We 
here employ the Born approximation, and the perturbed 
Green's function reads (7''^' = [z — iJ^'^*^'^ ±177]"^, where 
rj = is the magnitude of the self-energy. 

Before using the Green's function obtained above to 
calculate the longitudinal conductivity, the vertex cor- 
rection has to be considered. With the iterative method 
used in the case of the TI, we obtain the modified current 
operator j^, which is obtained by doubling the a in j^. 
The conductivity is then calculated via the Kubo formula 
as 

a.. = -4iTrj./jV = ^^ + ^^ (C2) 

with the r being the relaxation time, and no is the local 
electron density. The contribution of the first term of 
Eq. (jC2[) is the same as that of the Drude model, and 
we can see that the second term is the same as the con- 
ductivity of the surface state of the TI. The calculation 
of the spin accumulation can be done in a similar way. 
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which yields 



Ex, 



(C3) 



which is half the spin accumulation on the surface of the 
TI. 



Appendix D: Calculation by the Kubo formula 



where /3 = (A;bT)-\ i^n = W/3 (" = 0, ±2, ±4, • • • ) 
and w„ = mr/P in — ±1, ±3, ±5, • • • ). Subsequently, we 
make the analytic continuation ivm uj+i5 and take the 
limit limtj^o Qi{'-^)/i'-^ to obtain the response functions. 
The summation over n on the right hand side of Eq. (jPip 
is done by replacing the summation with the contour 
integral after multiplying that by the Fermi distribution 
function. With the aid of the residual theorem, we can 
represent the correlation function as 



In the previous sections, we adopted the Keldysh for- 
malism, which is convenient to deal with higher order 
perturbations. In this section, we perform the calculation 
of the spin polarization and the longitudinal conductivity 
of the TI via the Kubo formula^**. 

We perform the calculations by means of the Matsub- 
ara Green's function Q. To calculate the response func- 
tions of spin polarization (i — x,y,z), which are de- 
fined as ((Ti) — (j)iEx, we first compute the spin-current 
correlation function 



hn 



Az 



L2 / 27rie'3^~n"^'^^'^*^'°^^^*'^"'''^'^^'"*'^'*' ^^"^^ 



where we took the integral contour C as seen in Fig. [HI 
The integration can be further rewritten as 



(D3) 



r 



where f{z) = (e*^^ + 1) ^ is the Fermi distribution func- 
tion and —iAk{e) = Qk{^ + iS) — Qk{^ — 



Inn A 



77i/j3' 




contour<^ -InMp," 
-i V 



Re 



'branch cut 



FIG. 6: The contour of the integration. The contour is de- 
formed to avoid the branch cuts on the complex plane. 

By the analytic continuation ivm uj + i5 and Q{uj ± 
iS) = g^S^\ the correlation function Qi yields 
1 r de 



i{uj + iS) 



£2 / 2Tri 



{-fie) + fie + Lo)} 



(D4) 



We here consider the zero temperature limit where the 
Fermi distribution function becomes the step function. 



and hence {-f{e) + f{e + uj)}/u} -6{e) {oj 0). Then 
we obtain the response functions (jji as: 



which gives 4>x = 'Pz — and 



3'j:9k,t 



ev 



" "-imp ''imp 



(D5) 



(D6) 



These results are consistent with the results by the 
Keldysh formalism as expected. 

Similarly, we obtain the longitudinal conductivity by 
substituting cr,; with jx in Eq. (jPip . This leads to 



(^xx = -i^^ Tr jxgl,ejx9l,e^ , (D7) 



and finally we obtain 



TI'fT-impVJmp 



(D8) 



which is in accordance with the expression obtained by 
the Keldysh formalism. 
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